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ABSTRACT 

We investigate a stationary particle acceleration zone in the outer magnetosphere of 
a spinning neutron star. The charge depletion due to a global current causes a large 
electric field along the magnetic field lines. Migratory electrons and/or positrons are ac- 
celerated by this field to radiate curvature gamma-rays, some of which collide with the 
X-rays to materialize as pairs in the gap. As a result of this pair-production cascade, 
the replenished charges partially screen the electric field, which is self-consistently 
solved together with the distribution of particles and gamma-rays. If no current is 
injected at neither of the boundaries of the accelerator, the gap is located around 
the so-called null surface, where the local Goldreich- Julian charge density vanishes. 
However, we find that the gap position shifts outwards (or inwards) when particles 
are injected at the iimer (or outer) boundary. We apply the theory to the nine pulsars 
of which X-ray fields are known from observations. We show that the gap should be 
located near to or outside of the null surface for the Vela pulsar and PSR B1951-I-32, 
so that their expected GeV spectrum may be consistent with observations. We then 
demonstrate that the intrinsically large TeV flux from the outer gap of PSR B0540-69 
is absorbed by the magnetospheric infrared photons to be undetectable. We also point 
out that the electrodynamic structure and the resultant GeV emission properties of 
millisecond pulsars are similar to young pulsars. 

Key words: gamma-rays: observation ~ gamma-rays: theory - magnetic field - pul- 
sars: individual (B0540-69, B1055-52, B1509-58, B1951-I-32, Geminga, J0437-4715, 
J0822-4300, J1617-5055, Vela) 



1 INTRODUCTION 

The EGRET experiment on the Compton Gamma Ray Ob- 
servatory has detected pulsed signals from seven rotation- 
powered pulsars (for Crab, Nolan et al. 1993, Fierro et al. 
1998; for Vela, Kanbach et al. 1994, Fierro et al. 1998; for 
Geminga, Mayer-Hassel Wander et al. 1994, Fierro et al. 
1998; for PSR B1706-44, Thompson et al. 1996; for PSR 
B1046-58, Kaspi at al. 2000; for PSR B1055-52, Thomp- 
son et al. 1999; for PSR B1951-I-32, Ramanamurthy et al. 
1995). The modulation of the 7-ray light curves at GeV en- 
ergies testifies to the production of 7-ray radiation in the 
pulsar magnetospheres either at the polar cap (Harding, 
Tademaru, & Esposito 1978; Daugherty & Harding 1982, 
1996; Dermer & Sturner 1994; Sturner, Dermer, & Michel 
1995; Shibata, Miyazaki, & Takahara 1998), or at the vac- 
uum gaps in the outer magnetosphere (Cheng, Ho, & Ru- 
derman 1986a, b, hereafter CHR; Chiang & Romani 1992, 



1994; Romani and Yadigaroglu 1995; Romani 1996; Zhang 
& Cheng 1997, hereafter ZC97). Effective 7-ray production 
in a pulsar magnetosphere may be extended to the very high 
energy (VHE) region above 100 GeV as well; however, the 
predictions of fluxes by the current models of 7-ray pulsars 
are not sufficiently conclusive (e.g., Cheng 1994). Whether 
or not the spectra of 7-ray pulsars continue up to the VHE 
region is a question which remains one of the interesting 
issues of high-energy astrophysics. 

In the VHE region, positive detections of radiation at a 
high confidence level have been reported from the direction 
of the Crab, B1706-44, and Vela pulsars (Nel et al. 1993; 
Edwards et al. 1994; Yoshikoshi et al. 1997; see also Kifune 
1996 for a review), by virtue of the technique of imaging 
Cerenkov light from extensive air showers. However, as for 
pulsed TeV radiation, only the upper limits have been, as 
a rule, obtained from these pulsars (see the references cited 
just above) . If the VHE emission originates the pulsar mag- 
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netosphere, rather than the extended nebula, a significant 
fraction of them can be expected to show pulsation. There- 
fore, the lack of pulsed TeV emissions provides a severe con- 
straint on the modeling of particle acceleration zones in a 
pulsar magnetosphere. 

In fact, in CHR picture, the magnetosphere should be 
optically thick for pair-production in order to reduce the 
TeV flux to an unobserved level by absorption. This in turn 
requires very high luminosities of tertiary photons in the 
infrared energy range. However, the required fluxes are gen- 
erally orders of magnitude larger than the observed values 
(Usov 1994). We are therefore motivated by the need to con- 
trive an outer-gap model which produces less TeV emission 
with a moderate infrared luminosity. 

High-energy emission from a pulsar magnetosphere, in 
fact, crucially depends on the acceleration electric field, 
along the magnetic field lines. It was Hirotani and Shi- 
bata (1999a,b,c; hereafter Papers I, II, III), and Hirotani 
(2000a,b,c; hereafter Papers IV, V, VI) who first solved the 
spatial distribution of together with particle and 7-ray 
distribution functions. They considered a pair-production 
cascade in a neutron star magnetosphere by solving the 
Vlasov equations (see also Beskin et al. 1992). 

By this method, they explicitly solved the gap width 
along the field lines, E^^ , particle densities, and the 7-ray 
distribution functions. They further demonstrated that a 
stationary gap is formed around the null surface at which 
the local Goldreich-Julian charge density. 



PGJ 



2nc 



(1) 



vanishes (fig. |l|), where is the component of the mag- 
netic field along the rotation axis, SI refers to the angular 
frequency of the neutron star, and c is the speed of light. 
Equation (|^) is valid unless the gap is located close to the 
light cylinder, of which distance from the rotation axis is 
given by 



t^7LC = 3 X 10*5^2 ^cm, 



(2) 



where = r2/(100rad s~^). In this paper, we develop the 
method presented in Paper VI as follows: 

• We compute the curvature radius, Rc , at each point along 
the magnetic field line, rather than assuming Rc ~ O.S-tolc 
throughout the gap. 

• We partially take the effect of the unsaturated particles 
into account (eq.|38[). 

• We investigate the gap structure when particles fiow into 
the gap from the inner or the outer boundaries. 

• We compute the explicit spectra of TeV emission due 
to inverse Compton scatterings, rather than estimating only 
their upper limits. 

In the next two sections, we describe the physical pro- 
cesses of pair production cascade in the outer magnetosphere 
of a pulsar. We then apply the theory to individual pulsars 
and present the expected GeV and TeV spectra for various 
boundary conditions on the current densities in § ^ In the 
final section, we compare the results with previous works. 




Figure 1. Side view of a hypothetical outer-magnetospheric gap. 
On the null surface, the Goldreich-Julian charge density vanishes. 



2 BASIC EQUATIONS AND BOUNDARY 
CONDITIONS 

We first reduce the Poisson equation for the electric poten- 
tial into a one-dimensional form in 



2.1 



We then give the 
Bol tzma nn equations of the particles and the 7-rays in §§ 2.2 
and I2.3I a nd describe the Lorentz factors and the X-ray fields 
in §§ and 2.5. We finally impose suitable boundary con- 
ditions in 5 2.6. 



2.1 Reduction of the Poisson Equation 

In the Poisson equation for the electrostatic potential 
the special relativistic effects appear in the higher orders of 
{Qzu / c)'^ , where zu indicates the distance of the point from 
the rotational axis (Shibata 1995). As the first-order approx- 
imation, we neglect such terms. To simplify the geometry, 
we further introduce a rectilinear approximation. Unless the 
gap width, W , along the field lines becomes a good frac- 
tion of tULC, we can approximate the magnetic field lines 
as straight lines parallel to the s-axis (figure ^). Here, s is 
an outwardly increasing coordinate along a magnetic field 
line, while y designates the azimuth. Assuming that typical 
transfield thickness of the gap, D±, is greater than W , we 
can Fourier analyze in the transfield direction to rewrite the 
Poisson equation into (eq. [3] in Paper VI) 



$ : 



+ 47re iVj- 



N- 



PGJ 
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(3) 



where A''+ and refer to the positronic and electronic 
densities, respectively, e the magnitude of the charge on the 
electron, and s the length along the last-open fieldline. In 
what follows, we define the origin of s as the neutron star 
surface. 

It is convenient to non-dimensionalize the length scales 
by a typical Debey scale length c/cjp, where 



rric 2TTce 



1.875 X 10^J72^''^B5^'^^rad s"^ (4) 



Bs represents the magnetic field strength at the gap center 
in 10^ G unit. The dimensionless coordinate variable then 
becomes 
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Figure 2. Rectilinear approximation to the outer gap. The 
j/-axis, which designates the azimuth, is not depicted to avoid 
complication. This approximation is valid unless the gap width, 
W = S2 — si, becomes a good fraction of tolc- 



reaction forces cancel out each other to obtain the following 
continuity equations 



ds \ B J 



1 



de^ [rip+G+ + rip-G- 



(13) 



where G±{x,e^) are the distribution function of 7-ray pho- 
tons having momentum dznieCe^ along the poloidal field line. 
Since the electric field is assumed to be positive in the gap, 
e''''s (or e~'s) migrate outwards (or inwards). The pair pro- 
duction redistribution function for outwardly (or inwardly) 
propagating 7-rays are written as ryp+ (or rjp^) and are de- 
fined by 

'7p±(<^7) = (1 - Mc)— / rfex^^crp(e^,ex,Mc), (14) 
ujp ./ . dex 



C = {LUp/c)s = 6.25 X IQ-^Qi^/'^Bs^/^s. 



(5) 



By using such dimensionless quantities, we can rewrite the 
Poisson equation into 



p _ dip 



A, 



B. 



Ben 



(6) 



(7) 



where = {ujp/c)D± and v'(C) = e^{x) / {rricC''); the par- 
ticle densities are normalized in terms of a typical value of 
the Goldreich-Julian density as 



2-Kce N± {x) 
Q, B(x) 



(8) 



Bent = B(Ccnt) refers to the magnetic field strength at the 
gap center, ^ = ^cnt- In the non-relativistic limit, the local 
Goldreich-Julian charge density is given by 

B^C) _ 2cos6»cos(6'-ai)-sin6'sin(e'-ai) 

Bent \/l + 3 COS^ 6cnt 

where ai refers to the inclination angle of the magnetic mo- 
ment, 9 and Sent the colatitude angle of the position at which 
Pgj is measured and that at the gap center, respectively. Ne- 
glecting the general relativistic effect near to the star sur- 
face, we can implicitly solve in terms of s from 



i(6l - Qi)Y/l + 3cos2(6'-Qi) 



^LC .la sin6'LC sin^(0LC — Qi 



d9. 



(10) 



where 0* is the colatitude angle of the last-open fieldline at 
the star surface; it is given by 



sin2(6». 



sin^(6>Lc - Qj) 
■iz'Lc/sin6'LC 



(11) 



where r, refers to the neutron-star radius, and Slc the co- 
latitude angle of the intersection of the last-open fieldline 
and the hght cylinder. We can solve Slc for a given ai from 



2cot(SLC — Qi) sin6lLC + cos6'lc = 0. 



2.2 Particle Continuity Equation 



(12) 



Let us now consider the continuity equations of the parti- 
cles. It should be noted that almost all the particles migrate 
with the large, saturated Lorentz factor. We therefore as- 
sume here that both electrostatic and curvature-radiation- 



where ap is the pair-production cross section and cos~^/ic 
refers to the collision angle between outwardly (or in- 
wardly) propagating 7-rays and the X-ray photon with en- 
ergy mcC^tx', X-ray number density between dimensionless 
energies ex and e^+de^, is integrated over from the thresh- 
old energy eth = 2(1 — /ic)~^e7^ infinity. The explicit 
expression of ap is given by (Berestetskii et al. 1989) 



AcrT(l-«2) 



(3 _ In i±i: _ 2,;(2 - 1;') 
1 — V 



v{e^,, Ex, Mc) = a/ 1 ^ 



1 — e7ex 



(15) 



where ctt is the Thomson cross section, and ex = E^/mcC? 
refers to the dimensionless energy of an X-ray photon. 

To evaluate /Xc between the surface X-rays and the 7- 
rays, we must consider 7-ray's toroidal momenta due to an 
aberration of light. It should be noted that a 7-ray photon 
propagates in the instantaneous direction of the particle mo- 
tion at the time of the emission. The toroidal velocity of a 
particle at the gap center (rcnt,^cnt) becomes 



V4, = r-cnt sin Ocnt^ + kB^ 



(16) 



where k is a constant. On the right-hand side, the first term 
is due to corotation, while second one due to magnetic bend- 
ing. Since By arises in the gap, the corresponding drift ve- 
locity appears as the third term. Unless the gap is located 
close to the light cylinder, we can neglect the terms con- 
taining B0 as a first-order approximation. We thus have 
V4, ~ fcnt sin6'cntf2. Since the three-dimensional particle ve- 
locity is virtually c, the angle between the particle motion 
and the poloidal plane becomes 



-1 / Tent Sin6lcntf2 



1 / r-cnt Sm Ucnt 



p^h = sm 'I I = sm 'I I . (17) 

C / V Ci7LC / 

We evaluate (rcnt,^cnt) along a Newtonian dipole. Using (jj^h, 
we can compute /ic as 



= COS0ab COSfypol, 



(18) 



where Spoi is the collision angle between the surface X-rays 
and the curvature 7-rays, that is, the angle between the two 
vectors (rcnt,6'cnt) and (B'',B*) at the gap center. Therefore, 
the collision approaches head-on (or tail-on) for inwardly (or 
outwardly) propagating 7-rays, as the gap shifts towards the 
star. 

For magnetospheric, power-law X-rays, as assume 



© 2001 RAS, MNRAS 000, 



4 K. Hirotani and S. Shibata 



fic = COs(W/lI7Lc), 



(19) 



for both inwardly and outwardly propagating 7-rays. Aber- 
ration of light is not important for this component, because 
both the X-rays and the 7-rays are emitted nearly at the 
same place, provided that W <ti roLC holds. 

Let us introduce the following dimensionless 7-ray den- 
sities in the dimensionless energy interval between and 
f3i such that 

<?±'(C) = ^/ de,G±ix,e,). (20) 

In this paper, we set /3o = 10^; this implies that the lowest 
7-ray energy considered is 771^0^ — 51.1 MeV. To cover a 
wide range of 7-ray energies, we divide the 7-ray spectra 
into 9 energy bins and put /3i = 3.162 x 10^, (32 = 10^, /?3 = 
3.162 X 10^ 134 = 10", 05 = 3.162 x 10*, /Js = 5.623 x 10*, 
Pi = 10^ /3s = 1.778 X 10^ /39 = 3.162 x 10^ 

We can now rewrite equation ( [l3[ ) into the dimensionless 
form. 



(21) 



where rjp-j- are evaluated at the central energy in each bin as 

'^3^-l+P^\ ^22) 



A combination of equations ( pl| ) gives the current con- 
servation law. 



jtot = n+{^) + n_(^) = constant for ^. 



(23) 



When Jtot ~ 1.0, the current density per unit flux tube 
equals the Goldreich-Julian value, f2/(27r). 

2.3 Gamma- Ray Boltzmann Equations 

Let us next consider the 7-ray Boltzmann equations. The 7- 
rays are beamed in the direction of the magnetic field where 
it was emitted. Therefore, the propagation direction of each 
7-ray photon does not coincide with the local magnetic field 
where the 7-ray distribution is evaluated. However, to avoid 
complications, we simply assume that the outwardly (or in- 
wardly) propagating 7-rays dilate (or constrict) at the same 
rate with the magnetic field (Paper VII) . Then we can write 
down the 7-ray Boltzmann equations into the form. 



ds 



= rip±G±{x, 



where (e.g., Rybicki, Lightman 1979) 
ie^r 1 



-F 



hRc 
1 3 ftcr^ 



(!)• 



+ -r?cAf±(a;), (24) 
c 



(25) 



F{s) 



Rc 



K5 {t)dt; 



(26) 



(27) 



Rc is the curvature radius of the magnetic field lines and 
A's/s is the modified Bessel function of 5/3 order. The effect 
of the broad spectrum of curvature 7-rays is represented by 
the factor F{e^/e^) in equation (psj). In the polar coordinates 
{r,9), the explicit expression of _Rc is given by 



Rc 



-'I 



3/2 



d\ 



(28) 



where 

C = r-cos6»/ti7Lc, (29) 

C = rsinf/tULC. (30) 

Assuming a Newtonian dipole magnetic field with inclina- 
tion Qi, we can write the derivatives in equation(^) as 

d( cosoi — (3/2) sin&sin(6' — Oi) 



d^ sin ai -I- (3/2) cos ^ sin(6' — Qfi) ' 

d^C 1 (dC/rfC)/-sin6'cosai + (3/2) cos 6 sin a; 



Vi-|-(3/2) cose'sin(6' - Oi) 



where 

sin2(6l-Qi) _ sin2(e'. -Qi) 



(31) 
(32) 

(33) 



/ = (cos6'sincii — ^ sin^cosQi)^ -I- i cos(6 -I- Oi). (34) 

Integrating Jj^l ) in the energy intervals [Pi-i, /3i], we 
obtain 

± i9±i0 = ±i (InB) 9±' - Vp±9±' + r?^|^n±, (35) 
at; at, iJcnt 

where i = 1, 2, • ■ ■, m (m — 9) and 

^ ' du I Kr,{t)dt 




(36) 



2.4 Terminal Lorentz Factor 

The most effective assumption for particle motion in the 
gap arises from the fact that the velocity immediately sat- 
urates in the balance between the electric force and the ra- 
diation reaction force due to curvature radiation. Equating 
the electric force e\d<^/dx\ and the radiation reaction force, 
we obtain the saturated Lorentz factor at each point, 
1/4 



Tsat — 



3i?c^ d£ 
2e dx 



+ 1 



(37) 



When W is so small that a significant fraction of the par- 
ticles are unsaturated, equation ( p7| ) overestimates the par- 
ticle's Lorentz factors. To take account of such unsaturated 
motion of the particles, we compute F by 



1 

f " 
where 



(38) 



(39) 



represents the maximum attainable Lorentz factor; the gap 
is considered to exists in ^1 < ^ < ^2 and the potential origin 
is chosen such that ip{S,i) ~ 0. 
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2.5 X-Ray Field 

To execute the integration in equation (p^, we must specify 
the X-ray field illuminating the gap. X-ray field of a rotation- 
powered neutron star within the light cylinder can be at- 
tributed to the following three emission processes: 

(1) Photospheric emission from the whole surface of a cool- 
ing neutron star (Greenstein & Hartke 1983; Romani 1987; 
Shibanov et al. 1992; Pavlov et al. 1994; Zavlin et al. 1995). 

(2) Thermal emission from the neutron star's polar caps 
which are heated by the bombardment of relativistic parti- 
cles streaming back to the surface from the magnetosphere 
(Kundt & Schaaf 1993; Zavlin, Shibanov, & Pavlov 1995; 
Gil & Krawczyk 1996). 

(3) Non-thermal emission from relativistic particles accel- 
erated in the pulsar magnetosphere (Ochelkov & Usov 
1980a,b; El-Gowhari & Arponen 1972; Aschenbach & 
Brinkmann 1975; Hardee & Rose 1974; Daishido 1975). 

The spectrum of the first component is expected to be 
expressed with a modified blackbody. However, for simplic- 
ity, we approximate it in terms of a Plank function with 
temperature feTg, because the X-ray spectrum is occasion- 
ally fitted by a simple blackbody spectrum. We regard a 
blackbody component as the first one if its observed emitting 
area, A^, is comparable with the whole surface of a neutron 
star. At = 47rr*^, where r, denotes the neutron star radius. 
We take account of both the pulsed and the non-pulsed sur- 
face blackbody emission as this soft blackbody component. 
At distance r from the center of the star, the X-ray density 
between energies ex and tx -I- de^ is given by the Planck law, 
3 

/ As 



dex 



1 I meC" 
4^ 



Kinr^J exp(ex/A,) - 1' 



(40) 



where As is defined by 



(41) 

kTs refers to the soft blackbody temperature measured by 
a distant observer. Since the outer gap is located outside 
of the deep gravitational potential well of the neutron star, 
the photon energy there is essentially the same as what a 
distant observer measures. 

As for the second component, we regard a blackbody 
component as the heated polar cap emission if its observed 
emitting area, Ah, is much smaller than A^,. We approximate 
its spectrum by a Planck function. We take account of both 
the pulsed and the non-pulsed polar cap emission as this 
hard blackbody component. In the same manner as in the 
soft blackbody case, the spectrum at a radius r becomes 

3 



diVx 
dex 



1 

4^ 



ch 



(—) 

V 47rr2 J 



exp(ex/As) - 1' 



(42) 



where Ah is related with the hard blackbody temperature, 
fcTh, as 

AhEE^. (43) 

Unlike the first and the second components, a power-law 
component is usually dominated by a nebula emission. To 
get rid of the nebula emission, which illuminates the outer 
gap inefficiently, we adopt only the pulsed components of a 
power- law emission as the third component. We describe a 
magnetospheric component with the power law. 



dex 



pltx 



(en 



< ex < en 



(44) 



The photon index a is typically between —2 and —1 for 
a pulsed, power-law X-ray component in hard X-ray band 
(e.g., Saito 1998). We assume e^in = 0.1keV/511keV and 
emax ~ 100keV/511keV for homogeneous discussion. 



2.6 Boundary Conditions 

To solve the Vlasov equations (^, (^, (pl|), and (|35[), we 
must impose boundary conditions. The inner boundary is 
defined so that E« vanishes there. Therefore, we have 



S|i(Ci) =0. 



(45) 



It is noteworthy that condition ( [45|) is consistent with the 
stability condition at the plasma- vacuum interface (Krause- 
Polstorff & Michel 1985a, b) . At the inner (starward) bound- 
s-ry (C = Ci)i w6 impose (Paper VI) 



•^(Ci) = 0, 

5+^(6) =0 



1,2, 



(46) 
(47) 



where m — 9. Since the particles may be created inside of 
the outer gap (i.e., at ^ < ^i), positrons may fiow into the 
gap at ^ = ^1 as a part of the global current pattern in 
the magnetosphere. We thus denote the positronic current 
divided by ce per unit magnetic flux tube at ^ = as 

n+(6)=ii. (48) 

This yields, with the help of the charge-conservation law 
(eq. 11), 

n-(Ci)=itot-ji. (49) 

At the outer boundary (^ — ^2), we impose 

Si|(6)=0, (50) 

5.^(6) =0 (i = l,2,---,m), (51) 

n-(6)=i2. (52) 

The current (divided by ce) created in the gap per unit 
fiux tube can be expressed as 



Jgap = Jtot - Jl - J2. 



(53) 

We adopt jgap, ji, and j2 as the free parameters in this 
paper. 

We have totally 2m + 6 boundary conditions (U5[)-(|52) 
for 2m + 4 unknown functions $, n+, n_, g+ , g+ , 
■ ■ ■, (?+'", g-^ , g~^, ■ ■ •, fl-™'. Thus two extra boundary 
conditions must be compensated by making the positions of 
the boundaries ^1 and ^2 be free. The two free boundaries 
appear because = is imposed at both the boundaries 
and because jgap is externally imposed. In other words, the 
gap boundaries (^1 and ^2), and hence the position shifts if 
ji and/or j2 varies. 



3 TEV SPECTRA 

In TeV energies inverse Compton (IC) scatterings of infrared 
(IR) photons off relativistic electrons and positrons (P ~ 
10^) are the pro cess of 7-ray production. We briefiy discuss 
the IR field in 5 3.1, and describe the intrinsic TeV emission 
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from the gap in § 3.2 and the extrinsic absorption due to 



magnetospheric IR field in § B.J 



3.3 Absorption due to pair production 



3.1 Infrared photon field 

Without a careful consideration of the synchrotron beaming, 
it is difficult to estimate the specific intensity of the IR field. 
Therefore, we simply assume that the IR field are homoge- 
neous and isotropic within the radius zulc ■ This assumption 
comes from the following emission scenario in the magneto- 
sphere: The primary 7-rays emitted in the gap collide with 
X-ray and IR photons outside of the gap to materialize as 
primary pairs, which emit secondary photons in X-ray and 
soft 7-ray energies via synchrotron process. Some portions 
of the secondary photons collide with each other to mate- 
rialize as tertiary pairs in the magnetosphere. The tertiary 
pairs have much lower energies compared with the secondary 
pairs and emit copious tertiary IR photons via synchrotron 
process. On these grounds, we may expect that the tertiary 
IR field has lost the directional information of the primary 
7-rays, compared with the secondary X-ray field. 

We do not adopt a broken power-law to describe the 
soft photons from IR to X-ray energies. This is because the 
adopted assumptions are different each other. As described 
in § we assume that the secondary, power-law X-ray 
field is homogeneous only within the radius of rent around 
the gap and is highly beamed along the magnetic field at the 
gap cener so that the collision angles between the primary 
7-rays and the secondary X-rays in the gap are typically 
W/ui-Lc radian (eq. [^). On the other hand, we assume 
that the tertiary, power-law IR field is homogeneous and 
isotropic within the radius wlc- 

At distance r from the center of the star, we adopt the 
following power-law IR spectra: 



iViR 



d 
kpc 



(— ) 



(54) 



where emiricC refers to the IR photon energy, and tlR,min < 



We adopt eir,, 



10"" and eiR,„,i„ = 10" 



the results do not depend on these cut-off energies very 
much. 



3.2 Inverse Compton Scatterings 

When an electron or a positron is migrating with Lorentz 
factor r 1 in an isotropic photon field, it upscatters the 
soft photons to produce the following number spectrum of 
7-rays (Blumenthal & Gould 1970): 

d'^N 3 c dA'^iR deiR 

— — fJT 

dtd€j 4 deiR eiR 



2glng+(l + 2g)(l-g) + 



iQqfil-q) 
2(1 + Qq) 



(55) 



where Q = 4eiRr and q = e^/Q(r — e^); here, refers to the 
energy of the upscattered photons in nieC^ unit. Substitut- 
ing equation (|54|), integrating d^N/dtde^ over eiR, and mul- 
tiplying the 7-ray energy (e^meC^) and the electron number 
(A'e) in the gap, we obtain the fiux density of the upscattered 
photons as a function of e-y. 



As mentioned in 5 3.1, the IR field is assumed to be homo- 



geneous and isotropic within radius cc7lc. Outside of zulc, 
both the IR photon density and the collision angles decrease; 
therefore, we neglect the pair production at r > wlc for 
simplicity. The optical depth then becomes 



dNiE 
deiR 



■o"p(eiR, £7, Mc)deiR, 



(56) 



where L refers to the path length. For homogeneous dis- 
cussion, we assume L = 0.5ti7Lc in this paper. As the first 
order approximation, we adopt 90° as the collision angles in 
the isotropic IR field to substitute /jLc = in equation (p"" 
where tTp is defined by equation (n 



4 APPLICATION TO INDIVIDUAL PULSARS 

In this section, we apply the theory to the nine rotation- 
powered pulsars of which X-ray field at the outer gap can 
be deduced from observations. The Crab pulsar was inves- 
tigated in detail in Paper VII; therefore, we exclude this 
young pulsar in this paper. We first describe their X-ray 
and infrared fields in the next two subsections, and present 
the electric field distribution in 5 4.S and the resultant GeV 



and TeV emissions from individual pulsars in § 4.4 



4.1 Input X-ray Field 

We present the observed X-ray properties of individual pul- 
sars in order of spin-down luminosity, Eiot (table 1). We 
assume tmin = 0.1keV/511keV and e^ax = lOOkeV/SllkeV 
for homogeneous discussion. 

B0540-69 From ASCA observations in 2-10 keV band, 
its X-ray radiation is known to be well fitted by a power- 
law with a = —2.0. The unabsorbed luminosity in this 
energy range is 1.3 x lO'^^erg s~^, which leads to A^pi — 
9.0 X 10''d^(ro/ccLc)~^ (Saito 1998), where d is the distance 
to the pulsar in kpc. 

B1509-58 From ASCA observations in 2-10 keV band, 
its pulsed emission can be fitted by a power-law with 
a = —1.1. The unabsorbed flux in this energy range 
is 2.9 X 10""erg/cm^/s, which leads to A^pi = 9.3 x 
10"d2(ro/roLc)"^ 

J1617-5055 and J0822-4300 These two pulsars have 
resemble parameters such as f2 ~ 90rad s~^, n ~ 
-|^g30.6Q gj^3^ characteristic age ~ 8 x 10^ yr. From 

the ASCA observations of J1617-5055 in 3.5-10 keV band, 
its pulsed emission subtracted the background and the 
steady components can be fitted by a power-law with 
a = —1.6 (Torii et al. 1998). Adopting the distance to be 
3.3 kpc (Coswell et al. 1975), we can calculate its unab- 
sorbed fiux as 3.1 X 10~^^erg/cm'^/s, which yields A^pi = 
6.5 X 10^°d^{ro/zuLc)~^- On the other hand, the distance 
of J0822-4300 was estimated from VLA observations as 
d = 2.2 ± 0.3kpc (Reynoso et al. 1995). ROSAT observa- 
tions revealed that the soft X-ray emission of this pulsar is 
consistent with a single-temperature blackbody model with 
fcTs = 0.28 ± O.lOkeV and ~ 0.04^.(d/2.2)2 (Petre, 
Becker, & Winkler 1996). 
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Table 1. Input X-ray field 



pulsar 


distance 


f2 


logio M 




As /A, 




Ah/A, 


^pl 


—a 




kpc 


rad s 


Ig(G cm'^) 


eV 




cV 




cm ^ 




B0540-69 


49.4 


124.7 


31.00 










1014.15 


2.0 


B1509-58 


4.40 


41.7 


31.19 










1014.04 


1.1 


J1617-5055 


3.30 


90.6 


30.78 










1012.64 


1.6 


J0822-4300 


2.20 


83.4 


30.53 


280 


0.040 










Vela 


0.50 


61.3 


30.53 


150 


0.066 










B1951+32 


2.5 


159 


29.68 










1013.55 


1.6 


Geminga 


0.16 


26.5 


30.21 


48 


0.16 






105.00 


1.6 


B1055-52 


1.53 


31.9 


30.03 


68 


7.3 


320 


10-3.64 






J0437-4715 


0.180 


1092 


26.50 


22 


0.16 


95 


10-3.28 







Vela From ROSAT observations in 0.06-2.4 keV, the spec- 
trum of its point source (presumably the pulsar) emission 
is expressed by two components: Surface blackbody compo- 
nent with kTs = 150eV and = 0.066A,{d/0.5f , and a 
power-law component with a — —3.3 (Ogelman et al. 1993). 
However, the latter component does not show pulsations; 
therefore, we consider only the former component as the X- 
ray field illuminating the outer gap. 

B1951+32 From ROSAT observations in 0.1-2.4 keV, the 
spectrum of its point source (presumably the pulsar) emis- 
sion can be fitted by a single power-law component with a = 
— 1.6 and intrinsic luminosity of 2.3 x 10^^(ci/2.5)^erg s~^ 
(Safl-Harb & Ogelman 1995), which yields TVpi = 9.1 x 
10^^d'^(ro/TOLc)~^. The extension of this power-law is con- 
sistent with the upper limit of the pulsed component in 2-10 
keV energy band (Saito 1998). 

Geminga The X-ray spectrum consists of two compo- 
nents: the soft surface blackbody with kT^ = 50 eV and 
As = 0.22A* (d/0.16)^ and a hard power law with a — —1.6 
and A^pi = 3.9 x 10*' (Halpern & Wang 1997). A paral- 
lax distance of 160pc was estimated from HST observations 
(Caraveo et al. 1996). 

B1055-52 Combining ROSAT and ASCA data, 
Greiveldinger et al. (1996) reported that the X-ray spectrum 
consists of two components: a soft blackbody with fcTg = 68 
eV and As = 7.3^4* (d/1. 53)^ and a hard blackbody with 
fcTh = 320 eV and = 2.3 x 10"''4.(d/1.53)^ 
J0437-4715 Using ROSAT and EUVE data (Becker & 
Triimper 1993; Halpern et at. 1996), Zavlin and Pavlov 
(1998) demonstrated that both the spectra and the light 
curves of its soft X-ray radiation can originate from hot po- 
lar caps with a nonuniform temperature distribution and 
be modeled by a step-like functions having two different 
temperatures. The first component is the emission from 
heated polar-cap core with temperature fcTh — 10^'^^ K 
measured at the surface and with an area A^ = 5.3 x 
10"''yl.(d/0.180)^ The second one can be interpreted as a 
cooler rim around the polar cap on the neutron star sur- 
face with temperature kTs = 10^'^^ K and with an area 
As = 0.16j4* (d/0.180)^ . Considering the gravitational red- 
shift factor of 0.76, the best-fit temperatures observed at 
infinity become kTs = 22 eV and kT^ = 95 eV. From paral- 
lax measurements, its distance is reported to be 178 ± 26 pc 
(Sandhu et al. 1997). We adopt d = 180 pc as a representa- 
tive value. 



4.2 Input Infrared Field 

We next consider the homogeneous and isotropic infrared 
field in the magnetosphere. In addition to the references 
cited below, see also Thompson et al. (1999). 

B0540— 69 Its de-extincted optical and soft X-ray 
pulsed fiux densities can be interpolated as Fi, = 0.207i^~° '' 
Jy (Middleditch & Pennypacker 1985). This line extrapo- 
lates to 0.47 mjy at 640 MHz, which is consistent with an 
observed value 0.4 mJy (Manchester et al. 1993). We thus 
extrapolates the relation F^, — 0.207i^~'''^ to the infrared 
energies and obtain (Paper V) 

dNm/dem = 1.6 x lO^^d^em"^'^ (57) 

B1509— 58 The flux densities emitted from close to 
the neutron star in radio (Taylor, Manchester, & Lyne 1993), 
optical (Caraveo, Mereghetti, & Bignami 1994), soft X-ray 
(Seward et al. 1984), and hard X-ray (Kawai et al. 1993) 
bands can be fitted by a single power-law Fi, = 1.36u~^'^^ 
Jy. We thus adopt 

dNm/dem = 4.7 x lO^^d^eiR "^'^ (58) 

Vela The fiux densities emitted from close to the neu- 
tron star in radio (Taylor, Manchester, & Lyne 1993; Downs, 
Reichley, & Morris 1973), optical (Manchester et al. 1980) 
bands can be fitted by F^ = 1.71 x 10^1^"° ''^ Jy. We thus 
adopt 

dNm/dem = 1.8 x lO^d^em (59) 

B1951+32, The flux densities emitted from close to 
the neutron star in radio (Taylor, Manchester, & Lyne 1993) 
and soft X-ray (Safi-Harb & Ogelman, & Finley 1995) bands 
can be interpolated as Fi, — 32.8i/~'''*^ Jy. We thus adopt 

dNm/dem = 6.3 x lO^^d^em"^ '^ (60) 

Geminga The upper limit of the flux density in ra- 
dio band (Taylor, Manchester, & Lyne 1993) and the fiux 
density in optical band (Shearer et al. 1998) gives spectral 
index greater (or harder) than —0.69. Interpolating the in- 
frared fiux with these two frequencies, we obtain 

dNm/dem = 1.9 x lO^d^em"' ^ (61) 

which gives a conservative upper limit of infrared photon 
number density under the assumption of a single power-law 
interpolation. 

B1055— 52 The flux densities emitted from close to 
the neutron star in radio (Taylor, Manchester, & Lyne 1993) 
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and optical (Mignani et al. 1997) bands can be interpolated 



as = 4.2 X 10^1/"° " Jy. We thus adopt 
dNm/dem = 7.7 x lO^d^em"^ ■^ 



(62) 



J1617-5055, J0822-4300, and J0437-4715 There 
have been no available infrarad or optical observations for 
these three pulsars. We thus simply assume that a — —1.5 
for these four pulsars and that uF^ — 10^ Jy-Hz at 0.01 eV. 
We then obtain 

dNm/dem = 2.5 x lO^d'eiR"' ' (63) 
for J1617-5055, 

dNm/dem = 2.1 x lO^d^m"^ '^ (64) 
for J0822-4300, and 

dNm/dem = 3.6 x lO^d^iR"^ '^ (65) 
for J0437-4715. 

4.3 Electric Field Structure 

To examine the behavior of the solutions in a wide parameter 
space, we consider the following four cases: (^1,^2) = (0,0), 
(0.3,0), (0.6,0), and (0,0.3). In what follows, we denote 
them as cases 1, 2, 3, and 4, respectively. In case 1, par- 
ticles flow into the gap at neither the inner nor the outer 
boundary. In case 2 (or case 4), the positronic (or electronic) 
current flowing into the gap per unit flux tube at the outer 
(or inner) boundary is 30% of the typical Goldreich-Julian 
value, Q/2n. 

In case 1, for Qi — 45°, the solution for B1509-58 dis- 
appears if jgap > 0.0825 and that for B1055-52 disappears 
if igap > 0.0152. This is because the electric field distribu- 
tion forms a 'brim' (e.g., fig. 2 in Hirotani & Okamoto 1998) 
as jgap approaches a certain upper limit, typically several 
percent of Q/2n. On these grounds, we fix jgap = 0.01 for 
all the four cases. We consider B1509-58, B1055-52, and 
J0437-4715 to describe Eu distributions. 



4.3.1 B 1509-58, a young pulsar 

Since the Crab pulsar was investigated in Paper VII, we con- 
sider here B 1509-58 as a typical example of young pulsars. 
The results of -E||(s) for ai — 45° is presented in figure ^ 
The solid, dashed, dash-dotted, and dotted lines correspond 
to the cases 1, 2, 3, and 4, respectively. The abscissa desig- 
nates the distance along the last-open field line and covers 
the range from the neutron star surface (s = 0) to the posi- 
tion where the disance equals s = 1.2 x zui^c ~ 9.34 x 10® m. 

The solid line shows that the gap is located around the 
null surface, if no currents penetrate into the gap at neither 
of the boundaries. The gap slightly extends outwards, be- 
cause the pair production mean free path becomes small at 
the outer region of the gap, owing to the diluted X-ray field 
there. It is noteworthy that the potential drop becomes the 
largest when ji = and j2 = (i.e., case 1). 

The dashed and dash-dotted lines indicate three points 
we may notice. For one thing, the gap shifts outwards as 
ji increases. This conclusion is consistent with the results 
obtained analytically by Shibata and Hirotani (2000). What 
is more, the maximum value of decreases as the gap 



shifts outwards. This is because the decreased poj at larger 
distance reduces \dE^\/ds\ (eq.|^). One final point is that W 
increases as the gap shifts outwards. This is because the pair 
production mean free path increases, owing to the diluted 
X-ray field at large radii. The potential drop in the gap 
remains almost unchanged between cases 2 and 3, because 
the decrease of maximum and the increase of W cancel 
each other. 

On the other hand, when 72 increases, the gap shifts 
inwards. When the gap is located on the half way from the 
null surface to the star (case 4, the dotted line), the po- 
tential drop becomes only 3% of that in case 1. A physical 
interpretation on this point is given in 



5.2 



4. ,3. 2 PSR B1055-52, a middle-aged pulsar 

We next consider PSR B1055-52 as an example of middle- 
aged pulsars. The E\\ distribution for ai = 45° is presented 
in figure ^. The abscissa, ordinate, and the lines and the 
same as figure ^ 

Comparing with figure ^ we can understand that the 
gap is more extended for B 1055-52. This is because the X- 
ray field, which is supplied by the surface emission, becomes 
less dense for middle-aged pulsars. The weak X-ray field en- 
larges the pair-production mean free path, and hence the gap 
width. The strength of E\\ decreases compared with young 
pulsars, because local pcj is small due to their large tolc, 
or equivalently, due to their small f2. 

The position of the gap behaves in the same way as 
B1509-58. That is, it shifts towards the light cylinder (or 
the star surface) when ji (or J2) increases. 



4.3.3 JO437-4715, a millisecond pulsar 

We finally consider J0437-4715, a millisecond pulsar. The 
_E|I distribution for Qi = 45° is presented in figure |^. Even 
though the magnetic moment is small, the strength of E^^ is 
comparable with the young pulsar B1509-58. This is because 
the magnetic field in the outer magnetosphere is compara- 
ble with those of young pulsars, owing to its fast rotation 
{Q. = 1092rad s"^), which reduces iijlc to 2.74 x lO^m. In 
addition, W is much less than B1055-52, a middle-aged pul- 
sar and relatively close to B 1509-69, a young pulsar. This 
is because the shrunk light cylinder of millisecond pulsars 
results in a dense X-ray field illuminating the gap; the dense 
X-ray field in turn reduces the pair-production mean free 
path, and hence W. This result is consistent with the semi- 
analytical prediction for the millisecond pulsar B1821-24 
(Paper V), of which X-ray field is dominated by a power- 
law component. In short, the electrodynamic structures of 
the gap are similar between millisecond pulsars and young 
pulsars. 



4.4 Gamma-ray Spectra 

The spectra of outwardly (or inwardly) propagating 7-rays 
in GeV energies, are readily computed from g+^{(,2) (or 
(;_*(^i)). The TeV spectra, on the other hand, can be ob- 
tained by the method described in § ^. Except for B0540-69, 
the optical depth (eq. [Q) is much less than unity; there- 
fore, the effect of absorption due to TeV-eV collisions in the 
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31509-58 

inclination= 45 dec 




distance along field line [meter] 

Figure 3. Distribution of the acceleration field in the magneto- 
sphere of B1509— 58, a young pulsar. The inclination angle is 45°. 
The abscissa is the distance from the star surface along the last- 
open fieldline. The solid, dashed, dash-dotted, and dotted lines 
correspond to the cases 1, 2, 3, and 4, respectively (see text). 
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Figure 5. Distribution of the acceleration field for J0437-4715, 
a millisecond pulsar. The inclination angle is 45°. The abscissa, 
ordinate, and the lines are the same as figure figure 3. 



B1055-52 

inclination= 45 dec 




5x10" 10' 

distance along field line [meter] 

Figure 4. Distribution of the acceleration field for B1055— 52, 
a middle-aged pulsar. The inclination angle is 45°. The abscissa, 
ordinate, and the lines are the same as figure figure 3. 



B0540-69 

inclination= 45 deg 




frequency [Hz] 

Figure 6. Pulsed gamma-ray spectra from B0540-69 magneto- 
sphere when a\ = 45° . The thick and thin lines denote the spectra 
of outwardly and inwardly propagating 7-rays, respectively. The 
solid, dashed, dash-dotted, and dotted lines correspond to the 
cases 1, 2, 3, and 4, respectively (see text). 



magnetosphere, can be neglected for the other eight pulsars 
considered in this paper. 

We present the combined GeV-TeV spectra for B0540- 
69, B1509-58, J1617-5055, J0822-4300, Vela, B1951+32, 
Geminga, B1055-52, and J0437-4715. in figures mul- 
tiplying the same cross sectional area for the GeV and TeV 
emissions for each pulsar. Both the inchnation angle and 
D± are indicated at the top right corner. The outward (or 
inward) 7-ray flux is depicted by thick (or thin) lines. The 
observed pulsed fluxes and their upper limits are indicated 
by open circles (EGRET), open squares (Whipple), open tri- 
angles (Durham group). The upper limits of the stationary 
fluxes are obtained by CANGAROO group and are denoted 
by flUed circles. To compare the results with the same Qi, 
we present the case of ai = 45° for all the pulsars. For 
Vela, B195H-32, Geminga, and B1055-52, we also present 



the case of another ai to compare the expected spectra with 
EGRET observations. 



General Features 

We first consider the general features of the 7-ray spectra. 
T he G eV luminosity is the largest and hardest for case 1 
(§ 4.3), because both the potential drop and the maximum 

It should be noted that 



of E\\ is the largest, (see figs, 
the larger Oi is, the harder the GeV spectrum becomes; this 
can be understood if we compare figures ^ and for in- 
stance. It should be noted that the same conclusion was 
derived independently for a vacuum gap in Paper V (e.g., 
Ec in table 3), in which the 'gap closure' condition was con- 
sidered in stead of solving the stationary Vlasov equations 
(§ Physically speaking, this is because the distance of 
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Figure 7. Pulsed gamma-ray spectra from B1509-58 when «i 
45°. The Unes correspond to the same cases as in figure 7. 
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Figure 10. Pulsed gamma-ray spectra for the Vela pulsar when 
cti = 45° . The lines correspond to the same cases as in figure 7. 
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Figure 8. Pulsed gamma-ray spectra PSR J1617-30 when ai 
45° . The lines correspond to the same cases as in figure 7. 
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Figure 11. Pulsed gamma-ray spectra for the Vela pulsar 
when Oi = 75° . The dashed and dotted lines correspond to the 
same cases as in figure 7, while the dash-dot-dot-dotted lines to 

(igap,ii,i2) = (0.01,0.07, 0). 
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Figure 9. Pulsed gamma-ray spectra for the Vela pulsar when 
cti = 45° . The lines correspond to the same cases as in figure 7. 



the gap from the star decreases with increasing Qi for the 
same set of (jgap, 71,^2). At small radii, the strong magnetic 
field enlarges pgj, and hence \dE^^/ds\, which results in a 
large at the gap center, especially when the gap is nearly 
vacuum (as the solid lines indicate in figs. |l^ and [l^) . 

We can also understand that the TeV fluxes are gener- 
ally kept below the observational upper limits for appropri- 
ate GeV fluxes, or equivalently, appropriate cross sectional 
area, D±^. Therefore, the problem of the excessive TeV flux 
does not arise for reasonable IR densities of individual pul- 
sars. In the next subsection, we consider the expected spec- 
tra source by source. 

4-4-2 Comparison with EGRET observations 

First, we consider the Vela pulsar. It follows from the flg- 
ures |l^ and |ll| that the pulsed GeV spectrum obtained by 
EGRET (Kanbach et al. 1994; Fierro et al. 1998) is impos- 
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Figure 12. Pulsed gamma-ray spectra from B1951+32 when 
Oi = 45° . The lines correspond to the same cases as in figure 7. 
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14. Pulsed gamma-ray spectra for the Geminga pulsar 
= 45°. The lines correspond to the same cases as in 
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Figure 13. Pulsed gamma-ray spectra from B1951-f32 when 
o; = 75° . The lines correspond to the same cases as in figure 7. 

sible to be explained if the gap is located well inside of the 
null surface, as the dotted lines indicate. A large inclination 
angle and small current inflows are preferable to explain the 
observed fiux around 10^^ Hz. For example, we present the 
case for ji = 0.07 and j2 = as the (thick and thin) dash- 
dot-dot-dotted lines in figure For ai > 45° , expected TeV 
fluxes are well below the observed upper limit (Yoshikoshi 
et al. 1997). 

Secondly, we consider B1951-I-32. It follows again that 
the observed GeV flux around 10^^ Hz (Ramanamurthy et 
al. 1995) is impossible to be emitted if the gap is located 
well inside of the null surface, as the dotted lines indicate. 
Although the expected TeV flux is much less than the ob- 
served upper limits around lO^'^ Hz (Srinivanan et al. 1997) , 
it may be detectable around v ~ lO^'^ Hz with a future 
ground-based telescopes, if ai is as large as 75°. 

Thirdly, we consider the Geminga pulsar. For this 
middle- aged pulsar , only the dotted lines are depicted in fig- 
ures ^ and ^ because there exist no solutions for cases 1, 
2, and 3. That is, for ai = 45° (or 60°), solutions exist only 
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Figure 15. Pulsed gamma-ray spectra for the Geminga pulsar 
when a; = 45°. The lines correspond to the same cases as in 
figure 7. 



when j2 > 0.073 (or j2 > 0.017), provided that jgap = 0.01 
and ji = hold. This is because its X-ray field is so weak 
that the gap extends towards the light cylinder and vanishes 
dE\^/ds at ^ = ^2 (i.e., forms a 'brim') at the critical value 
of j2 (= 0.073 for Ui = 45°, say). Since the solutions exist in 
a limited region of the parameter space, the observed GeV 
spectrum above lO'^'* Hz (Mayer-Hasselwander et al. 1994) 
are impossible to be explained by the present theory. The 
TeV flux is expected to be much less than the observed up- 
per limits (for pulsed upper limits, see Akerlof et al. 1993; for 
steady flux and/or its upper limits, see Weekes and Helmken 
1977, Kaul et al. 1985, Vishwanath et al. 1993, Bowden et 
al. 1993). 

Fourthly, we consider B1055-52. For this middle-aged 
pulsar, solutions exist in a wide parameter range, unlike the 
Geminga pulsar. However, we obtain too soft GeV spectra to 
flt the observations (Thompson et al. 1999). The difficulties 
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Figure 16. Pulsed gamma-ray spectra from B1055— 52 when Figure 19. Pair production optical depth for a TeV photon to 

Oi = 45° . The lines correspond to the same cases as in figure 7. be absorbed in the homogeneous, isotropic IR field in the mag- 

netosphere of pulsar B0540-69. The abscissa designates the 7-ray 

frequency. 
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Figure 17. Pulsed gamma-ray spectra from B1055— 52 when 
Qi = 75° . The lines correspond to the same cases as in figure 7. 
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Figure 18. Pulsed gamma-ray spectra for J0437-4715 when 
Qi = 45° . The lines correspond to the same cases as in figure 7. 



5.1 



Finally, we consider the other pulsars of which GeV 
fluxes were not detected by EGRET. For B0540-69, the ba- 
sic properties are common with the Crab pulsar (Paper Vll) 
as follows: 

(1) W <C ta^lc holds because of its large X-ray density. 

(2) Intrinsic TeV flux is comparable or even somewhat 
greater than the GeV flux but is significantly absorbed by its 
dense, magnetospheric IR field. For example, in case 1, the 
intrinsic TeV flux before absorption attains 1.7 x 10^'' Jy Hz 
at 1.4 X 10^^ Hz, which is about two times greater than the 
maximum of the GeV flux. In figure ^ we present the ab- 
sorption optical depth (eq. Q). The resultant, absorbed 
TeV spectra are depited in figure ^ together with the GeV 
spectra. 

(3) The curvature spectra are relatively hard (fig. be- 
cause its strong magnetic field at the gap increases _E|| , al- 
though its W is small. 

For B1509-58, its relatively less dense IR fleld (com- 
pared with Crab and B0540-69) cannot absorb enough TeV 
fiux in the magnetosphere. As a result, the expected TeV 
fiux becomes comparable with the observational upper limit 
while the GeV flux is well below the observational upper 
limits (flg. 0). As flgures ^ and |l| indicate, the TeV fiux 
does not increase so rapidly as GeV one with increasing Qi. 
We thus conclude that ai should be comparable with or 
greater than 45° so that the GeV and TeV fluxes may be 
consistent with the observations. (Note that we can adjust 
the cross sectional area to decrease or increase the GeV and 
TeV fluxes at the same rate.) 

In the case of J1617-5055, its GeV spectrum (flg. ^) is 
softer than Crab (Paper Vll) and B0540-69, because its X- 
ray field is less dense compared with these two young pulsars. 
Its TeV fiux is small because of its small IR fleld assumed 

(§0)- 

The GeV spectrum of J0822-4300 (flg. g) is softer than 
that of Vela, although the X-ray fleld more or less resembles 
each other. This is because the former's high-temperature 
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surface emission provides a denser X-ray field, which shrinks 
the pair-production mean free path, and hence W to reduce 
, compared with Vela. 
For J0437-4715, the GeV flux is large (fig.Jl|), although 
we assume a small cross sectional area, D± = {1.5W)^. 
This large Gev flux comes from the fact that millisecond 
pulsars have similar electrodynamic structure with young 
pulsars. Compared with the EGRET upper limit of pulsed 
emissions (Thompson 2000), we can expect that this milli- 
second pulsar is an attractive target for the next-generation 
GeV telescopes. 



5 DISCUSSION 

In summary, we have developed a one-dimensional model for 
an outer-gap accelerator in the magnetosphere of a rotating 
neutron star. When the positrons (or electrons) penetrate 
into the gap from the inner (or outer) boundary, the gap 
shifts outwards (or inwards). Applying the theory to Vela 
and B1951-I-32, we flnd that the GeV spectrum softens sig- 
nificantly and cannot explain the EGRET observations, if 
the gap is located well inside of the null surface. We can 
therefore conclude the gap should be located near to or out- 
side of the null surface for these two pulsars. We also ap- 
ply the theory to B0540-69, B1509-58, J1617-5055, J0822- 
4300, and J0437-4715 and find that the TeV fluxes are un- 
detectable by current ground-based telescopes for moderate 
cross sectional areas of the gap for all the seven pulsars men- 
tioned just above, if we estimate their magnetospheric IR 
fleld by interpolating radio and optical pulsed fluxes. For 
B0540-69, its intrinsically large TeV flux is absorbed by 
its dense IR field to become undetectable; this conclusion 
is analogous with what was obtained in Paper VII for the 
Crab pulsar. 



5.1 Difficulties in tlie Application to Middle-aged 
Pulsars 

We briefly discuss the difficulties in the application to 
middle-aged pulsars. As figures |l4|-|l6| indicate, the fluxes 
from Geminga and B 1055-52 above 10 Hz are quite insuf- 
flcient compared with EGRET observations. There are three 
reasons for this problem. First, special relativistic effects be- 
comes important for the extented gaps of middle-aged pul- 
sars, of which X-ray fleld is weak. For example, in equa- 
tion (^), we neglected the terms containing higher orders 
in (Qzxj/c)^. Secondly, the global geometry of the magnetic 
fleld lines becomes important as as the gap extends. That 
is, the rectilinear approximation (§ 2.1) cannot be applied 



to these two middle-aged pulsars. Thirdly, their large 
which is a good fraction of the trans-field electric field, pre- 
vents particles from corotation (eq.[^). On these grounds, 
applications of the present theory to middle-aged pulsars do 
not provide good agreement with observations. 

5.2 Gap Closure Condition 



In 



4.3.1 



, we pointed out that the potential drop decreases 
significantly as the gap shifts inwards from the null sur- 
face. In this subsection, we interpret this behavior from the 
gap closure condition investigated in Papers IV and V. In 



a stationary gap, the pair production optical depth, W/\p, 
equals the ratio A'^(jgap/jtot), where Ap and A^.^ refer to the 
pair production mean free path and the number of 7-rays 
emitted by a single particle, respectively. We thus obtain the 
following gap closure condition: 



W - 



J gap 
y • 

Jtot 



(66) 



This condition is automatically satisfied by the stationar- 
ity of the Boltzmann equations. For example, consider the 
case when the external current is ten times greater than the 
created one in the gap. In this case, the expectation value 
for the 7-rays that are emitted from a single particle in the 
gap, should be one tenth so that a stationary pair production 
cascade may be maintained, because not only the produced 
particles but also the externally injected particles emit 7- 
rays in the same way. 

As the gap shifts inwards, the X-ray density increases 
to reduce Ap. At the same time, the ratio jgap/jtot decreases 
as _72 increases. As a result, W decreases very rapidly with 
increasing j2. Owing to the rapidly decreasing W , the inte- 
grated potential drop significantly decreases, although the 
local PGJ, and hence \dE^^/ds\ increases inwards. 

If the gap shifts outwards, on the other hand, the in- 
crease of Ap and the decrease of jgap/jtot partially cancel 
each other. As a result, the potential drop decreases only 
slightly with increasing ji . 



5.3 Comparison with ZC97 

It is worth comparing the present method with ZC97, who 
onsidered that the gap size is limited so that the emitted 
7-rays may have energies just above the pair production 
threshold in collisions with the surface X-rays that are due 
to the bombardment of produced particles in the gap. Their 
method qualitatively agrees with the current gap closure 
condition (eq. |6^), provided jgap/jtot ~ 1 (see § 5.8 in Pa- 
per V); case 1 (ji = j2 = 0) corresponds to this case. How- 
ever, as j2 increases, the factor (jgap/jtot) in equation ( |66| ) 
becomes small to reduce W . In other words, the closure con- 
dition adopted by ZC97 would become inconsistent with sta- 
tionary Vlasov equations (§ 2), if we were to overextraporate 
it to the case when j2 ~ 1. In the same manner, we cannot 
overextraporate the gap-closure condition of ZC97 to the 
case when ji increases and hence the gap shifts towards the 
light cylinder (as cases 2 and 3), because jgap/jtot decreases 
substantially. 



5.4 Inverse Compton Scatterings in Young Pulsar 
Magnetosphere 

Although the dense IR field of young pulsars (B0540-69, 
B1509-58) results in a relatively large intrinsic TeV emission 
compared with other pulsars, the intrinsic TeV flux is not 
more than the GeV flux, provided ai > 45°. Therefore, the 
radiation reaction force is caused primarily by the curvature 
process rather than the inverse Compton scatterings, as long 
as j2 is less than 0.3. However, if j2 exceeds 0.3, we cannot 
in general rule out the possibility of the case when the radia- 
tion reaction is caused by IC scatterings. This is because the 
dense X-ray field will suppress the particle Lorentz factors 
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(Paper II) below lO'^, and because such less-energetic par- 
ticles scatter copious IR photons into lower 7-ray energies 
with large cross sections (~ ctt). As a result, the GeV emis- 
sions may be mainly due to IC scattering in this case. More- 
over, since the magnetic field becomes very strong close to 
the star, not only the synchro-curvature process (ZC97), but 
also a magnetic pair production becomes important. There 
is room for further investigation of the case when an 'outer' 
gap is located close to the polar cap. 
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